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It was recently demonstrated [P. B. Weichman, Phys. Rev. Lett. 91, 143908 (2003)] that the 
scattered electric field from highly conducting targets following a rapidly terminated electromagnetic 
pulse displays a universal t~^^'^ power law divergence at early time. It is now shown that for strongly 
permeable targets, ^ic/ Hb ^ 1, where is the background magnetic permeability, the early time 
regime separates into two distinct power law regimes, with the early-early time t~^^^ behavior 
crossing over to t~^^'^ at late-early time, reflecting a spectrum of magnetic surface modes. The 
latter is confirmed by data from ferrous targets where jj.c/ fJ-b = 0(10^), and for which the early- 
early time regime is invisibly narrow. 



PACS numbers: 03.50.De, 41.20.-q, 41.20.Jb 

Remote characterization of buried targets is a key goal 
in many environmental geophysical applications, such as 
landmine and unexploded ordnance (UXO) remediation 
. A tool of choice is the time-domain electromagnetic 
(TDEM) method, in which an inductive coil transmits 
low frequency [typically 0(100 Hz)] EM pulses into the 
ground. Following each pulse [terminated rapidly on a 
ramp timescale — ©(lO^/is)], the voltage V{t) induced 
by the scattered field is detected by a receiver coil. Stan- 
dard TDEM sensors are capable of resolving signals from 
very small (of order 1 gram) metal targets |3], and are 
therefore well suited to UXO detection. 

Low frequency yields increased sensitivity to conduct- 
ing targets, and increased exploration depth (below 5 
m) but leads to nearly complete loss in spatial resolu- 
tion. Lacking direct spatial imaging capability (enabling 
straightforward identification), one is reduced to seek- 
ing such information indirectly via a careful analysis of 
the full time dependence of V{t). As described in Refs. 
0,0, this signal is affected by both intrinsic (target size, 
shape, geometry, and other physical characteristics) and 
extrinsic (relative target-sensor position and orientation, 
transmitter and receiver coil geometries, pulse waveform, 
etc.) properties, and the key to discrimination is the ex- 
traction of the former from the "background" of the lat- 
ter. The aim of this letter is to further develop such 
formalism for the early time part of V{t) |3]. 

Conductor electrodynamics are essentially diffusive, 
and the basic time scale Tc — L^/Dc is determined by 
the target diameter Lc, and diffusion constant Dc — 
c^/47r/Xc<Tc (in Gaussian units) depending on the (rela- 
tive) target permeability fic and conductivity CTc- Ferrous 
targets (e.g., steel) are typically modelled as paramag- 
netic targets with very large permeability = O(IO^), 
and (MKS) conductivity ctc = 0(10"^ S/m) . Thus, even 
for targets as small as = 1 cm, one finds decay times 
Tc — 0(10^ ms), much larger than typical pulse peri- 
ods. Larger UXO-like targets, which are usually ferrous, 
have even larger Tc. It will often be the case, therefore, 
that the full measured range of V{t) will lie in the early 
time regime, t <C Tc. On the other hand, for a nonmag- 



netic (e.g., aluminum) target of similar size one obtains 
Tc = 0(1 ms), and the measured signal will cover a much 
broader dynamical range. 




FIG. 1: (COLOR) Schematic diagram of measurement and 
early time dynamics. Left: prior to pulse termination, trans- 
mitter coil currents generate a magnetic field in the neigh- 
borhood of the target. Right: just after pulse termination, 
the target interior has not had time to adjust to the absence 
of the transmitted field, and screening surface currents are 
generated to enforce the correct EM boundary conditions. 

The early time response is governed by the dynamics of 
the screening currents, induced in response to the rapid 
quenching of the transmitted magnetic field immediately 
following pulse termination (see Fig.^). The initial dif- 
fusion of these currents inward from the target surface 
generates a t~^^^ power law in V{t), with coefficient re- 
flecting the surface properties of the target Q. However, 
underlying this power law is the assumption that t/Tc is 
the smallest parameter in the problem. It will now be 
shown that the background-target permeability contrast 
Mb/Mc can generate a new small parameter that greatly 
limits the its range of validity. Specifically, an extended 
calculation is described that divides the early time regime 
t < Tg, where Tg <^ Tc is the point at which bulk effects 
first begin to enter (see further below) , into an early-early 
time regime, < t ^ T^ag, where the V{t) ^ t~^/^ re- 
mains valid, and a late-early time regime Tmag ^ t ^ Tg, 
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where a new power law V{t) ~ t obtains. In the 
neighborhood of the magnetic crossover time, 

Tmag^Tcifib/Hcf, (1) 

an interpolation between the two power laws occurs (that 
is visible only if Tmag <C Tg), for which the full functional 
form is provided. For ferrous targets T,nag/Te = 0(10^''), 
and the early-early time regime is essentially invisible, 
and only the behavior is seen, consistent with mea- 
sured data: see Fig. |21 
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FIG. 2: (COLOR) Geonics EM61-3D data (black curves) from 
nonmagnetic (aluminum) and ferrous (steel) targets, together 
with predicted exact solutions 5] (blue curves). Asymptotic 
early-time power laws (red curves) are seen, with t~^^'^ for 
aluminum, t~^^'^ for steel. The former shows a crossover to 
multi-exponential decay behavior (magenta curve — the first 
60 slowest decaying modes), consistent with = £9(10~^s), 
Tc — 0(10^^8). The latter covers essentially the entire ob- 
served range, consistent with Te — O(lO^^s), Tc = 0(1 s). 

At low frequencies the dielectric function in the ground 
and in the target is dominated by the its imaginary part, 
e = A-Kia/uj, where (t(x) is the dc conductivity, and the 
Maxwell equations may be reduced to a single equation 
for the vector potential, 

„ / 1 \ 47r(7 „ 47r 

Vx -VxA +— atA=— js, (2) 
\M / c c 

with magnetic induction B = V x A, field H = B//i, and 
gauge chosen so that the electric field is E = — (l/c)9tA. 
The transmitter loop generates the source current den- 
sity j 5 (x, t). The conductivity and permeability are sepa- 
rated into background [crb(x), /iti(x)] and conducting tar- 
get [crc(x), /^c(x)] components, where ac, fJ-c vanish out- 
side the target volume Vc- High conductivity contrast, 
CTb/cTc ^ I, is required, but ^ib/lJi-c is arbitrary. 

Equation |(2Jl is a vector diffusion equation, with con- 
trast Db/Dc = 0(10^-10^). The "background commu- 
nication time" between instrument and target separated 
by distance i? is r;, = P? jD^. For representative val- 
ues (Tb = O.I S/m, /ib — 1, = 3 m, one obtains 



Tb = C'(IO~^s), instantaneous even on the scale of Tr- 
Thus, the electrodynamics of the background may be 
treated as quasistatic: outside the target one may drop 
the dtA. term in (j^J. Further, following pulse termina- 
tion {t > 0) one has js = 0, hence V x H = 0, and 
one may derive H = — V^> from a scalar potential $, 
which in turn satisfies V • (/ibV^) = 0. The external 
field is therefore entirely determined by the internal field 
through the boundary condition Q n- [/XbHb — /-tcHc] = 0, 
where Hb, He are the fields just outside and inside 
the target surface, and n is the local surface normal. 
Thus, one obtains a Neumann-type boundary condition, 
— n • V<i> = (/ic//ib)n • He, and hence the formal solution 

$(x)= / dVgjv(x,r)^n.H,(r), (3) 

where dVc is the target surface, and is the Neumann 
Green function satisfying — /i^^V • (iib'^gN) — i5(x— x') 
with boundary condition n • V^at = 0. 

The initial surface screening current appears in the 
transverse magnetic boundary condition 4, 5j: K(r) ~ 
(c/47r)n X [Hb(r) - Hc(r)]. As indicated in Fig.Q] He is 
the same as that prior to pulse termination, while Hb is 
given by (0): 

K = -^nx (V<f>-He)t=o, (4) 

which, via determines K entirely in terms of He. 

To investigate the evolution of K for i > 0, wc take 
advantage of the rapid variation of the fields near the 
surface with the local vertical coordinate z. Thus, z- 
derivatives dominate ||2Jl, and to leading order in the 
small parameter e{t) = \J DctjL^ 'g^, one obtains 

-a|A-L + a* A = 0, z< 0, (5) 

with initial condition E(0) = -(l/c)5tA(0) = ct^^K(5(z). 
Here Z = z/\/TTc, A-^ = A — n(n • A) is the tangential 
part of A, and /i, a, D are treated as constants on either 
side of the boundary. One may choose A A — for i = 0, 
and it follows immediately that n- AA = 0: A A is purely 
transverse. 

One requires now the boundary condition for Q at 
Z = 0". For t > let A = A(0)-t-AA(t) and $ = $(0)-t- 
A$(t), where one treats | AA|/|A(0)|, | A$|/|$(0)| = 
0{e) and |(92AA|/|AA| = 0{l/e). Keeping only lead- 
ing terms 6], continuity of n x H (the surface current 
sheet now has finite thickness) and Q imply that 



y^^K = -9zAA^|2=o- +Aic/D^nx VA$|z^o+, 

(6) 

in which A<i> is given by AA(Z = ) via ©. In fact, 

1 /2 

estimating dz = 0[Dc /e{t)Lc], gw — 0{1/Lc), and nx 
V — 0{1/Lc), the ratio of the second term on the right 
hand side to the first is 0[€{t)Lir/ fib]- Their relative order 
is therefore time- dependent Q. In particular, at early- 
early time e{t) <C iJib/ tJ^c (i-e., t <C Tmag) one may drop the 
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second term. Lack of time dependence in K then leads 
to the simple homogeneous Neumann boundary condition 
9zE = for the electric field. This is the limit in which 
the early time analysis in Ref. Q was carried out and the 
behavior of V{t) derived. 

Thus, for ^c/ = one gains nothing by in- 

cluding the A$ term: it is of the same order as other 
C'[e(i)] terms previously dropped from © 6]. However, 
if Hc/fJ'b S> 1 this term dominates for Tmag < t < Te, 
which now comprises a large fraction of the early time 
regime. The remainder of this paper is concerned with 
extending the theory into this regime. 

With the A$ term, © is nonlocal, coupling AA(Z — 
0~) over the surface. In order to decouple © we diago- 
nalize it by seeking transverse vector eigenfunctions a„, 
n = 1, 2, 3, . . ., satisfying 



K„a„(r) = -ficVDcii x V7/;„(r) 



(7) 



where k„ is the corresponding eigenvalue, and r is the 
surface coordinate. The scalar is defined by 



(8) 



where, to condense the notation, we define the 
Neumann Green function operator via (f>{r) — 
Jgy d'^r'gN{r,r')(j>{r'). Substituting into one ob- 
tains a scalar eigenvalue problem. 



(9) 



where the generalized surface Laplacian is £A</'(r) = n • 
V X [/ic VlJ^n X V0(r)] . On a sphere one obtains Ca = 
— Ur\/ DrL^^T? , where Lc is now the radius, and L = 



— ix X V is the angular momentum operator 

As the OLn are derived from the scalar il)n via Q , they 
do not form a complete set of 2D vector fields. The miss- 
ing fields consist of the kernel of Q (the space of func- 
tions K„ = 0), spanned by all fields (3 with vanishing 
normal magnetic field, n • V x /3 = 0. It follows that 
/3 = —V^cj) = n X n X V(t> for some other scalar Let 
/3ri be chosen as eigenfunctions of the generalized 
transverse Laplace equation. 



A„(/)„(r) = -n • V X 
/3„(r) - -V^0„, 



(10) 



with a new set of eigenvalues A„ . 

The set {Q:„,/3n,} forms a complete basis, and we per- 
form the magnetic surface mode expansion 



oo 

AA(x,i) = ^[ylW(Z,i)a(r)+A(2)(Z,t)/3(r 



n=l 



oo 

K(r) = ^[xWa(r)-|-i^(2)/3(r) 



(11) 



From 0, £^/.t(7^ri-VxAA = X)™^" Vn, and from dTUIl, 



n- V X (/x^ ^Dc ^^^AA) = ^„ XnAn'(^n (and similarly for 



K), hence appropriate orthogonality relations for (j>n 

may be used to determine An^ , Kn ^ . 

Substituting l(TT|) into and complete separation 
of the surface modes is achieved: 

{dt ~ dl)A[^ = 0, Z < 0, (12) 

with initial and boundary conditions 

(az + «;„<5.i)A«U=o- = -^i^^- (13) 
The solutions are An^ [Z, t) = —Kn^H{Z, t; KnSn), where 



H(Z,t;K) 
iJ(Z,t;0) 



erfc 



erfc 



— /4tp-^ /4t 



iZlerfc 



/ 2Kt-Z \ 

\ Vit ) 



(14) 



where erfc(a;) is the complementary error function. 

Fields external to the target are obtained by extending 
q:„,/3„ into the exterior space. First, let 



'0n(x) 
0n(x) 



— / d2^'5iv(x,r')-£AV'«(r') 



/ dV'gz5(x,r') 

JdV^ 



(15) 



in which (7£)(x, x') is the Dirichlet Green function sat- 
isfying — cTjJ^^V • (cTfeVgu) — (5(x — x') and vanishing on 
dVc, and x is no longer restricted to the surface. These 
definitions guaranteed continuity at the boundary. The 
vector eigenfunctions are now extended by solving 

^bV X a„(x) = -V?A„(x), V • {(JbOLn) = 0, (16) 

while imposing continuity of q:;|; at the surface, and 

/3„(x) = -V0„(x). (17) 

Since V x /3„ = it does not contribute to the external 
magnetic field, and hence to any inductive measurement. 
The external field is now simply 

oo 

AA(x,i) = [4'H0-,i)a„(x) +42)(0-,t)/3„(x) 

(18) 

with time-dependence given by the internal field on the 
boundary, hence governed by the function 



H{Q,t-K) 



1 



1 - e"" *crfc 



(19) 



/ii [1„ l(7rK2^)l/2+0(K2^)] ^ k2^«1 
(7rAt2t)-l/2+0[(K2i)-3/2]} ^2^» 1. 

The time derivative, entering the electric field and volt- 
age, displays the promised t^^/^ and t^'^/'^ power laws 
in the two limits. If '0„ varies on scale L„ < Lc, one 
may estimate from (O k„ = C'[rniag^(ic/in)], hence 
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FIG. 3; (COLOR) Comparison between exact and early time results for spheres in a homogeneous background. Plotted is 
Hi{t = t/Tc) for 1 < Z < 5 (curves for higher values of I lie lower on the plots). Left: iJ,c/lJ.b ~ 1. Blue (dashed) lines show 
the asymptotic early-early-time form I19II : red (dashed) lines show the full early time solution; black (solid) lines show exact 
results; and green (dash-dotted) lines show H21|l truncated at the first 3 terms (an approximation that might emerge from a 
late time perturbative approach 0]). The two early time curves, though differing in detail, exhibit roughly the same accuracy, 
with interval of validity shrinking as Lim/Lc = 1/^^, as predicted. In realistic applications, the transmitted field in the target 
region will be fairly uniform, and one may expect Z = 1 to dominate, with small corrections from higher I. The union of early 
and late time approximations then provide a rather accurate description of the full signal. Center: fic/l-ib ~ 5. The full early 
time form now clearly exhibits greatly extended accuracy over the early-early time power law. Right: /ic//.i6 = 100. The full 
early-time form is indistinguishable from the late-early time form II19II . the early-early time interval being invisibly narrow. 



crossover point Ti„ag,n = ^/k-I = Ohnaglin/ic)^]- One 
expects Ln = Lc only for the fundamental mode, hence 
([T)l actually represents an upper bound on the spectrum 
of crossover times Tmag,™- 

The exact solution for a homogeneous sphere in a ho- 
mogeneous background serves to clarify all of the above 
concepts. The operators and £a now commute and 
may be simultaneously diagonalized using spherical har- 
monics: CgYijn = Yim/{l + l)Lc, CAYlm = -[ficVD^l{l + 
l)/Ll]Yi„,,. Thus, using Vim = (f>hn = L^'^Yim one ob- 

— 1/2 

tains = iTmag (hcncc mode length scale L;„j — 
Lc/l), and Xim = / L'^IJ-cVT^c- The vector functions 

are Q;m(x) = -inb\/{l + I) /1{L{/ X^+^)'Xim, f3lm{^) = 
^i^{l + l)/lV X [iL[/x'+^)Ximi where X/^ = [1(1 + 
l)]^i/^LF;„ are the vector harmonics |^. The exact so- 
lution associated with the magnetic modes otim is gov- 
erned by the usual bulk exponentially decaying mode 
expansion |3, H, from which one obtains ^[^(i) ~ 



Ki2V^c[Hi{0) - Hi{t/T,)], where 

J/ (On) - Ji + l(On)ji-l(C/n) Ci„ 

where j;(x) are the spherical Bessel functions, and the 
scaled decay rates are given by the roots of 

= ^CW/-l(On) +lU-^ JliCln). (21) 

In Fig. ^ Hi{t/Tc) is compared to the early time pre- 
diction Hi{t/Tc) w Hi{0) - {4Tc)-^/^H{0-,t;Kim), with 
Hi{0) = {p.b/2^ic)/[l + {1 + l)^-lb/|J■c], for various I, Hc/fJ'b- 
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